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1. Introduction 

Let T > be a finite time horizon and T , P) be a stochastic basis supporting a d-dimensional 
Brownian motion W . We assume that the filtration F = (Tt)t<T generated by W satisfies the 
usual assumptions and that Tt = ^. 

Let (AT, Y, Z) be the solution of the decoupled Brownian Forward-Backward SDE 

Xt = Xo+ [ b{Xs)ds+ f <j{Xs)dWs (1.1) 
Jo Jo 

Yt = g{T,Xr)+J^ ls<rf{Xs,Ys,Z,)ds- Z.dWs ,ie[0,T], (1.2) 

where t is the first exit time of {t,Xt)t<T from a cylindrical domain D = [0,T) x O for some 
open piecewise smooth connected set O C M'', and b, a, f and g satisfy the usual Lipschitz 
continuity assumption. 

This kind of systems appears in many applications. In particular, it is well known that it is 
related to the solution of the semi-linear Cauchy Dirichlet problem 

— Cu — f{-,u, Dua) ^ on D , u = 5 on dpD , (1-3) 
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where C is the (paraboHc) Dynkin operator associated to X, i.e. for tjj € C^'^ 
C^j := dtiJ + {b, DV) + ^Tr [aD^] , a := aa* , 

and dpD ([0,T) x dO) U ({T} x O) is the parabohc boundary of D. More precisely, if the 
solution u of (1.3) is smooth enough, then Y = u{-, X) and Z = Dua{-, X). Thus, in the regular 
frame, solving (1.2) is essentially equivalent to solving (1.3). 

In this paper, we study an Euler scheme type approximation of (1.1)-(1.2) similar to the one 
introduced in [4, 29], see also [2, 3, 24]. We first consider the Euler scheme approximation X 
of X on some grid tt := {ti = ih, i < n} with modulus h := T/n, n £ N* . The exit time r is 
approximated by the first discrete exit time f of {ti, Xti)tien from D. Then, the backward Euler 
scheme of (F, Z) is defined for i = n — 1,...,0 as 



Yu := E [Yu^, I .Ft.] + lu<f h f{Xu , Yu , K 



Zt. /i-iE [y*.^, (m.+i 



with the terminal condition Yt = g(f,Xf) ■ Here, g is a suitable extension of the boundary 
condition on the whole space [0,T] x R'^. 

The main purpose of this paper is to provide bounds for the (square of the) discrete time 
approximation error up to a stopping time 6' < T P — a.s. defined as 



Err(/?,)^ maxE 


sup lt<e\Yt-Yu\'' 


+ E 


[L 




tG[t,,ii + i] 







\z. 



Z^(t)\\'^dt 



(1.4) 



where (j){t) := sup{.s G tt : s < t}. 

We are interested in two important cases: 6 = T and 9 = t Af. The quantity Err(/i)T coincides 
with the usual strong approximation error computed up to T. The term FjTT{h)rAf should be 
more considered as a weak approximation error, since the length of the random time interval 
[0,r A f] cannot be controlled sharply in pratice. It essentially provides a bound for Yq — Yq, 
or equivalently in terms of (1.3), u(0,Xo) — Yq. Let us mention that a precise analysis of the 
weak error has been carried out by Gobet and Labart in [14] in the uniformly elliptic case with 
O = W^. 

As in [4], [23] and [29], who also considered the limit case O = M.'^ (i.e. t = T), the approximation 
error can be naturally related to the error due to the approximation of X by X^ and the 
regularity of the solution (Y, Z) of (1.2) through the quantities: 



n{YYs. 



maxE 



sup 



\Yt - Yt. 



and 'R{Z)li^2 



where 



Zgds I Tt, 



for i < n . 



Z^{t)\\'^dt 



(1.5) 
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In the case / 0, y is a martingale and It; is the best approximation of Yt on the time 
interval by an JFt^ -measurable random variable. In this case, Doob's inequalities imply 



that E 



for 



universal constant c > 0. 

Moreover, the definition (1.5) impHes that is the best approximation in L^([0,T] x V.,dt(g) 
dF) of Z by a process which is constant on each time interval [ti,ti^i). Thus, TZ{Z)'!^2 < 

This justifies why TZ(Y)g2 and TZ{Z)'^2 should play a crucial role in the convergence rate of 
Err(/i) to as ^ 0. 

Bounds for similar quantities have previously been studied in [4, 29] in the case O = and 
in [2, 24] in the case of reflected BSDEs. All these articles use a Malliavin calculus approach to 
derive a particular representation of Z . Due to the exit time, these techniques fail in our setting. 
We propose a different approach that relies on mixed analytic/probabilistic arguments. Namely, 
we first adapt some barrier techniques from the PDE literature, see e.g. Chapter 14 in [11] and 
Section 6.2 below, to provide a bound for the modulus of continuity of u on the boundary, 
and then some stochastic flows and martingale arguments to obtain an interior control on this 
modulus. Under the standing assumptions of Section 2, it allows to derive that TZ{YYg2 + 
TZ{Z)'!^2 = Oih) and that u is l/2-H61der in time and Lipschitz continuous in space. 

To derive our flnal error bound on Err(/i)e, we additionally have to take into consideration the 
error coming from the approximation of r by f. We show that E [|t — f |] = 0{h^~^) for all 
e > 0. Combined with the previous controls on TZ{YYg2 and TZ{Z)'^2, this allows us to show 
that Err(/i)T = 0{h^~^). Exploiting an additional control on a weaker form of error on r — f, 
we also derive that Err(/i)T-AT = 0{h^~'^). As a matter of facts, the global error is driven by 
the approximation error of the exit time which propagates backward thanks to the Lipschitz 
continuity of u. 

Importantly, we do not assume speciflc non degeneracies of the diffusion coefficient but only a 
uniform non characteristic boundary condition and uniform ellipticity close to the corners, recall 
that O is piecewise smooth. Using the transformation proposed in [19], these results could be 
extended to drivers with quadratic growth (for a bounded boundary condition g). Also, without 
major difficulties, our results could be extended to time dependent domains and coefficients (6, 
cr and /) under natural assumptions on the time regularity. We restrict here to the homogeneous 
cylindrical case for simplicity. 

We note that the numerical implementation of the above scheme requires the approximation 
of the involved conditional expectations. It can be performed by non-parametric regression 
techniques, see e.g. [15] and [22], or a quantization approach, see e.g. [1] and [7, 8]. In both cases, 
the additional error is analyzed in the above papers and can be extended to our framework. We 
note that the Malliavin approach of [4] cannot be directly applied here due to the presence of 
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the exit time. Concerning a direct computable algorithm, we mention the work of Milstein and 
Tretyakov [25] who use a simple random walk approximation of the Brownian motion. However, 
their results require strong smoothness assumptions on the solution of (1.3) as well as a uniform 
ellipticity condition. 

The rest of the paper is organized as follows. We start with some notations and assumptions in 
Section 2. Our main results are presented in Section 3. In Section 4, we provide a first bound 
on the error: it involves the error due to the discrete time approximation of r by f and the 
regularity of the solution {Y, Z) of (1.2). The discrete approximation of t is specifically studied 
in Section 5. Eventually, Section 6 is devoted to the analysis of the regularity of (1.3) and (1.2) 
under our current assumptions. 

2. Notations and assumptions 

Any element a; G M**, d > 1, will be identified to a line vector with i-th component and 
Euclidean norm The scalar product on R'^ is denoted by {x,y). The open ball of center x and 
radius r is denoted by B{x, r), B{x, r) is its closure. Given a non-empty set A C M'', we similarly 
denote by B{A,r) and B{A,r) the sets {a; e M'' : d{x,A) < r} and {x e M.*^ : d{x,A) < r} 
where d{x, A) stands for the EucHdean distance of x to A. For a (m x d)-dimensional matrix M, 
we denote M* its transpose and we write M S M'^ if m = d. For a smooth function f{t, x), Df 
and D^f stand for its gradient (as a line vector) and Hessian matrix with respect to its second 
component. If it depends on some extra components, we denote by dtf(t, x, y, z), dxf{t, x, y, z), 
etc... its partial gradients. 

2.1. Euler scheme approximation of BSDEs 

From now on, we assume that the coefficients of (1.1)-(1.2) satisfy: 

(HL): There is a constant L > such that for all (t, x, y, z, t' , x', y', z') e ([0, T] x R'' x R x R'')^: 

\\{b,a,g,f){t,x,y,z)^ib,a,g,f){t',x',y',z')\\ < L \\{t, x,y, z) - {t' ,x' ,y\ z')\\ , 
\\{b,<j,gj){t,x,y,z)\\ < L (1 + ||(a;, y, z)||) . 

Under this assumption, it is well known, see e.g. [27, 28], that we have existence and uniqueness 
of a solution {X, Y, Z) in 5^ x 5^ x , where we denote by the set of real valued adapted 
continuous processes ^ satisfying ||^||52 := E [sup^^y ^ < oo , and by Ti? the set of 
progressively measurable R'^-valued processes C for which HCIIw^ := E[Jq'" jCipdiJ^ < oo . 

As usual, we shall approximate the solution of (1.1) by its Euler scheme X associated to a grid 

t: ■= {t'i = ih , i < n} , h :~ T/n , n G N* , 
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defined by 

Xt - ^0+/ b{X^(,))ds+ f cT{X^(,))dWs ,t>Q, (2.1) 

Jo JQ 

where we recall that (/)(s) := argmax{ti, i < n : ti < s} for s > . 

Regarding the approximation of (1.2), we adapt the approach of [29] and [4]. First, we approx- 
imate the exit time r by the first exit time of the Euler Scheme {t,Xt)teTT from D on the grid 
tt: 

f := inf{t e TT : ^ O} A T . 

Remark 2.1. Note that one could also approximate r by f inf{t e [0,T] : Xt 4- '^l ^ ^) 
the first exit time of the "continuous version" of the Euler scheme (t, -?t)te[o.T] , as it is done 
for linear problems, i.e. / is independent of (F, Z), see e.g. [13]. However, in the case where O 
is not a half-space, this requires additional local approximations of the boundary by tangent 
hyperplanes and will not allow to improve our strong approximation error, compare Corollaire 
2.3.2. in [12] with Theorem 3.1 below. 

Then, we define the discrete time process (F, on tt by 

Ft. := E[Ft.^J .Ft.] +1*.<^ /i/(Xt,,Ft,,Z*J , (2.2) 
Zt. := /i-^E [Ft.^, {Wu^^ - Wt^) I Tt^ ,i<n, (2.3) 

with the terminal condition 

Ft = gif,X^). (2.4) 
Observe that Yt^\ti>T = g{T,Xf)'^ti>T and that Zt^'\.ti>T = 0. 

One easily checks that (Fj. ,^4.) S for all i < n under (HL). It then follows from the 
martingale representation theorem that we can find Z £ Ti? such that 

Ft.+, - E [Yt^^, I Tu] = ZsdWs for all i<n. (2.5) 

This allows us to consider a continuous time extension of F in defined on [0, T] by 

Yt = g{f,Xr)+ [ W /(X^(.),F^(.),%s))ds- / ZsdWs. (2.6) 



Remark 2.2. Observe that Z = on ]r, T] and Z = on ]t, T]. For later use, also notice that 
the Ito isometry and (2.5) imply 

rti+i 



E 



Zsds I Ft, 



ti 



, i<n . (2.7) 
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2.2. Assumptions on C?, cr and g 

Our main result holds under some additional assumptions on O, cr and g. Without loss of 
generality, we can specify them in terms of the constant L which appears in (HL). 

We first assume that the domain O is a finite intersection of smooth domains with compact 
boundaries: 

(Dl): We have O := C\7=i where meW and is a domain of for each I < £ < m. 
Moreover, has a compact boundary, sup{||a;|| : x e dO^} < L, for each 1 < £ < m. 

It follows from Appendix 14.6 in [11] that there is a function d which coincides with the algebraic 
distance to dO, in particular O := {a; e R'' : d{x) > 0} , and is outside of a neighborhood 
B{C,L^^) of the set of corners 



We also assume that the domain satisfies a uniform exterior sphere condition as well as a uniform 
truncated interior cone condition: 

(D2): For all x G dO, there is y{x) G O", r{x) e [L-\L] and S{x) G 5(0, 1) such that 

B{y{x),r{x))nO = {x} and {x' e B{x, L'^) : {x' - x, S{x)) > {1 - L-^)\\x' - x\\} C O . 

In view of (Dl), these last assumptions are actually automatically satisfied outside a neighbor- 
hood of the set of corners, see e.g. Appendix 14.6 in [11]. 

In order to ensure that the associated first boundary value problem is well posed in the (uncon- 
strained) viscosity sense, we shall also assume that 



satisfies a non-characteristic boundary condition outside the set of corners C and a uniform 
eUipticity condition on a neighborhood of C: 

(C): We have 



m 



C := Pi 90^ n do'' . 



a := era 



mi{n{x)a{x)n{x)* 



X e dO\B{C,L-^ 



)} > L ^ where n{x) := Dd{x) 



and 



inf{^a(2;)^* 



^ e 95(0, 1) , X G O n B(C, L-^)} > L 



-1 



In particular, it guarantees that the process X is non-adherent to the boundary. 
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Observe that n coincides with the inner normal unit on dO outside the set of corners. By abuse 
of notations, we write n{x) for Dd{x), whenever this quantity is well defined, even if a; ^ dO. 

Importantly, we do not assume that cr is non degenerate in the whole domain. 

We finally assume that g is smooth enough: 

(Hg): g e Ci'2([0,T] x R'^) and \\dtg\\ + \\Dg\\ + \\D'^g\\ < L on [Q,T] x R'^ . 

Clearly, this smoothness assumption could be imposed only on a neighborhood of dO. Since it 
is compact and Y depends on g only on dO, we can always construct a suitable extension of g 
on M.'^ which satisfies the above condition. Actually, one could only assume that g is Lipschitz 
in (t, x) and has a Lipschitz continuous derivative in x. With this slightly weaker condition, all 
our arguments would go through after possibly replacing 5 by a sequence of regularized versions 
and then passing to the Hmit, see Section 6.4 for similar kind of arguments. 



3. Main results 

We first provide a general control on the quantities in (1.4) in terms of TZ{Y)'^2, TZ{Z)'!^2 and 
|r — f|. Let us mention that this type of result is now rather standard when O = M.'^, see e.g. 
[4], and requires only the Lipschitz continuity assumptions of (HL). 

Proposition 3.1. Assume that (HL) and (Hg) hold. Then, there exist Cl > and a positive 
random variable S^l satisfying E [(^l)^] < for all p > 2 such that 



Err(/i)^ < CL[h + niY)^2+n{Z);^2+E 



(3.1) 



and 



Err(ft)2^^ < Err(/i)2^^- < Cl (/^ + 7^(r)5. + 7^(Z)^0 + E E [a|r - f| | .F,^^^] 



Cl 



lr<rE 



n 2 



\Zs\\ds I Tf 



(3.2) 



where r+ is the next time after t in the grid tt; t+ := inf{t e tt : t < t} 



The proof will be provided in Section 4 below. Note that we shall control Err(/i)^^^ through 
the sHghtly stronger term Err(/i)^^^^, see (3.2). This will allow us to work with stopping times 
with values in the grid tt which will be technically easier, see Remark 4.2 below. 

In order to provide a convergence rate for Err(/i)y and Err(/i)^^^^, it remains to control the 
quantities TZ{Y)^2, 'R{Z)^2 and the terms involving the difference between t and f. 
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The error due to the approximation of r by f is controlled by the following estimate that extends 
to the non uniformly elliptic case previous results obtained in [12], see its Corollaire 2.3.2. The 
proof of this Theorem is provided in Section 5 below. 

Theorem 3.1. Assume that b and a satisfy (HL) and that (Dl) and (C) hold. Then, for 
e G (0,1) and each positive random variable satisfying E[(^)p] < C£ for all p > I, there is 
C£ > such that 



E 



E |t - f| I J^r 



< Cfh 



'E 1, 1-e 



r+ArJ 

In particular, for each s G (0, 1/2), there is C£ > such that 

E[|T-f|] < ^ 

In [12], the last bound is derived under a uniform ellipticity condition on a and cannot be ex- 
ploited in our setting, recall that we only assume (C). Up to the e term, it can not be improved. 
Indeed, in the special case of a uniformly eUiptic diffusion in a smooth bounded domain, it has 
been shown in [16] that E [r — f] = Ch^ + o{hi) for some C > 0, see Theorem 2.3 of this 
reference. 



Our next result concerns the regularity of (Y, Z) and is an extension to our framework of similar 
results obtained in [23], [4], [3] and [2] in different contexts. 



Theorem 3.2. Let the conditions (HL), (Dl), (D2), (C) and (Hg) hold. The 
Moreover, for all stopping times 6,-d satisfying 6 < ■& < T V — a.s., one has 



E 



and 



E 



sup \Y,^Yg 

?<i5<)? 



\ZsFds I 



2p 



< E[eL\d-9\ I J^g] , p-1,2. 



(3.3) 



(3.4) 



(3.5) 



where ^£ is a positive random variable which satisfies E < oo , for all q > 1. 

In addition, the unique continuous viscosity solution u of (1.3), in the class of continuous solu- 
tions with polynomial growth, is uniformly 1/2-IIdlder continuous in time and Lipschitz contin- 
uous in space, i.e. 



\u{t,x)-u{t',x')\<CL(\t~t'\i + \\x-x'\\^ forall{t,x) and{t',x') 



e D 



(3.6) 
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The proof is provided in Section 6 below. The bound (3.5) can be interpreted as a weak bound 
on the gradient, whenever it is well defined, of the viscosity solution of (1.3). It implies that Y is 
l/2-H61der continuous in norm. This result is rather standard under our Lipschitz continuity 
assumption in the case where O — M."^, i.e. t — T, but seems to be new in our context and under 
our assumptions. The bound TZ{Z)1^2 < Cl h can be seen as a weak regularity result on this 
gradient. It would be straightforward if one could show that Duu is uniformly l/2-H61der in 
time and Lipschitz in space, which is not true in general. 

Combining the above estimates, we finally obtain our main result which provides an upper 
bound for the convergence rate of Err(/i)^ (and thus for Err(/i)^^^) and Err(/i) 



JT- 



Theorem 3.3. Let the conditions (HL), (Dl), (D2), (C) and (Hg) hold. Then, for each 
£ G (0, |), there is Cf^ > such that 



Err(/i)2^^^ < Cl h^-" and Err(/i)| < C£ /i^"" . 

This extends the results of [2, 3, 29] who obtained similar bounds in different contexts. 

Remark 3.1. When r can be exactly simulated, we can replace f by t in the scheme (2.2)- 
(2.3). In this case, the two last terms in the right hand-sides of (3.1) and (3.2) cancel and we 
retrieve the convergence rate of the case O = R'', see e.g. [4]. 



Remark 3.2. Note that the Lipschitz continuity assumption with respect to the x variable on 
g and / is only used to control at the right order the error term coming from the approximation 
of X by X in g and /. If one is only interested in the convergence of Err(/i)T this assumption 
can be weakened. Indeed, if we only assume that 

(HL'i): b, a satisfy (HL), sup{|/(-, 2)|, {y, z) e M x R**} and g have polynomial growth, and 
f{x,-) is uniformly Lipschitz continuous, uniformly in x € R'', 

a weak version of (3.1) can still be established up to an obvious modification of the proof of 
Proposition 4.2 below. Namely, there exists C > and a positive random variable ^ satisfying 
E [{S.Y] < Cl for all p > 2 for which 



Err(/i) 



< C \ h+E 



CE 



\giT,X^)-gif,X^) 



-7^(Z)^2 +E 



\fiX,,Y,.Z, 



e|T-r|- 
~/(X^(,),n,Z,)|^ds 



'ds 



The terms 



Ytj)(s)\'^ds and TZ{Z)'!^2 are easily seen to go with h, see e.g. the proof 
of Proposition 2.1 in [.3] for details. As for the other terms in the first fine, it suffices to appeal to 
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Theorem 3.1 which impUes that E [^\t — r|] ^ and that r — > r in probabihty under (Dl) and 



(C). Note that the last assertion impUes that E 1^^^ J^J ||Zs|p(is 
in probabihty. Hence, under the additional continuity assumption 

(HL'2): g and f{-,y, z) are continuous, uniformly in {y, z) e K x M^, 

we deduce that the two last terms in the second Hne go to as well. 



and Xr - 



4. Euler scheme approximation error: Proof of Proposition 
3.1 

In this section, we provide the proof of Proposition 3.1. We first recall some standard controls 
on X, (y, Z) and X which holds under (HL). 

From now on, denotes a generic constant whose value may change from line to line but which 
depends only on Xq, L and some extra parameter 77 (we simply write Cl if it depends only on Xq 
and L). Similarly, denotes a generic non-negative random variable such that E — ^1'^ 

for all p > 1 (we simply write £,l if it does not depend on the extra parameter rj) . 



Proposition 4.1. Let (HL) hold. Fix p > 2. Let d he a stopping time with values in [0,T]. 
Then 



E 



sup 

t6[i?,T] 



\Ytr 



\Zt\\ dt 



and 



E 



sup {\\Xt\ 
te[i9,r] 



\Xtr) I 



< ci{i + \\x,r) 



Moreover, 



maxE 

z<n 



sup {\\Xt - Xt^P + \\Xt ~ Xu\\P) 



te[ti,ti 



E 



sup \\Xt-Xt\\P 
te[o,T] 



siip\\Xt- X^(^t)\\ >r 

t<T 



<CLr-^ h , r > 



and, if 9 is a stopping time with values in [0, T] such that ■d<9<i3 + h¥ — a.s., then 
E[\\Xe~X4'' + \\Xe-X4P \ J^^] < eLh^ ■ 
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Remark 4.1. For later use, observe that the Lipschitz continuity assumptions (HL) ensure 
that 



E 



sup \\Ytr- 

te[i?,T] 



\Ztfdt I 



< oo for all p > 2 . 



In order to avoid the repetition of similar arguments depending whether we consider Err(/i)0 
with 9 = T or 9 = T+ At, we first state an abstract version of Proposition 3.1 for some stopping 
time 9 with values in tt. 

Proposition 4.2. Assume that b,a and f satisfy (HL). Then, for all stopping time 9 with 
values in n, we have 



Errihfe < Cl { h + E [\Ye - Ye\^] + 7^(F)52 + 7^(Z)^2 + E 



r{tVT)Ad 
/ (a + lf<rl|^.!|')rf.' 



Let us first make the following Remark which will be of important use below. 



Remark 4.2. Let d < 9 P — a.s. be two stopping times with values in tt and H be some 
adapted process in S^. Then, recalling that ti^i — ti ^ h, it follows from (2.7) and Jensen's 
inequality that 



E 



H4,{s)\\Z4,{s)\\'^ds 





E 











Zudu I Tt^ 



< 



ML 



ti 

HtM-' / \\Z4^duds 
'ti Jti 



ds 



< E 



H^is)\\Zsrds 



By definition of Z, see (1.5), the same inequality holds with {Z, Z) or {Z — Z,Z — Z) in place 
of (Z, Z). This remark will allow us to control \\Z — Z^\\ through \\Z — Z\\ and \\Z — Z^\\, see 
(4.3) below, which is a key argument in the proof of Proposition 4.2. Observe that the above 
inequality does not apply if 1) and 9 do not take values in tt. This explains why it is easier to 
work with r+ instead of r, i.e. work on Err(/i)^^^^ instead of Err(/i)^/^^. 



Proof of Proposition 4.2. We adapt the arguments used in the proof of Theorem 3.1 in [4] to 
our setting. By applying Ito's Lemma to (Y — F)^ on [t A 9, ti^i A 9] for t G [ti, ti+i] and i < n, 
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we first deduce from (1.2) and (2.6) that 



E 



\Yt^0 — ^tAeP 



\Zs - ZsW ds 



tf\e 



{Ys - Ys) (i.<r/(e,) - i.<r/(e^(.))) ds 



where the martingale terms cancel thanks to Proposition 4.1 and Remark 4.1, and where := 
{X,Y,Z) and 6 := {X,Y,Z). Using the inequality 2ab < + b^, we then deduce that, for 
a > to be chosen later on, 

rti+iAd 

\Ys-Ys\^ds 



A 



< 



^[\Yu^^Ae-%^^As\^] +a: 



tA9 



ti + iA0 

tAB 
ti + iAB 

tAd 



i.<. (/(e.)-/(e^(,)))'ds + 



tAB 



ifmfds 



l-T<S<f 



Recall from Remark 2.2 that Z = on ]t,T]. Since Yt = g{T, Xr) on {t > t}, we then deduce 
from (HL) and Proposition 4.1 that 

pti+iAe 

E[\Yt^^,^g-%^,^e\'']+aE 



Cl a-^E 



Cl a"'E 



Cl a-^E 



h\Yt^Ae-%Ae? + 



tA0 
tAO 



\Ys~Z\^ds 



ti + 1 AS 



tAO 
ti + iAd 



h + \\Zs - Zms) 



Za,i_^)\\ ds 



tA0 



{£,l'^tAt<s<t\Jt + lr<s<r|l^s|l^)c^S 



(4.1) 



It then follows from Gronwall's Lemma that 

E^YtAe-YtAG?] < il + C"lh)E[\Yu^,^e-Yu^,Ae\' 



+ {Cl a-^ +C2 h)E 
+ {Cl +C2 h)E 



ti+iAB 



h lYt^AO-Yt^^el' + / |n-y^(,)|'ds 



tAB 



tA0 



h + \\Zs — Z^^s)\\ + \\Z(f,{s) — Z^(_s}\\ j ds 



+ {Cl +Cl h)E 



ti + iAd 



{£,L'i-TAf<s<TVf + If <s<r II ■^s|P)c's 



tAO 



(4.2) 
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Plugging (4.2) in (4.1) applied with t = ti, using Remark 4.2, taking a > large enough, 
depending on the constants Cl, and h small leads to 



< {l + CLh)V.[\Yu^,;,e~%+M''] 



tiAe 



ds 



{£.L'i-TAf<s<TWf + lr<s<T||^s|P)c?S 



tiAe 



This implies that 



— ZJ\ ds 









Z<p{s) - Z^^s)\\'^ds 


+ E 


/ \\Zs - Z^i^s)\\'^ds 






Jo 



< Cl {E[\Ye-Ye\']+h + n{Y)l,+n{Z);,,) 

+ ClE iL\fhe-TAe\+ u<,^r\\Zsfds 

Jo 

We conclude the proof by using Remark 4.2 again to obtain 
E ^ < Cl 

which imphes the required result, by the definition of Err(/i)g in (1.4) 
The above result implies the first estimate of Proposition 3.1. 

Proof of (3.1) of Proposition 3.1. It suffices to apply Proposition 4.2 for 6 = T and observe 
that the Lipschitz continuity of g imphes that 









Zs - ZsW'^ds 


+ E 


/ \\Zs - Z^i^s-^lfds 






Jo 



(4.3) 



□ 



E[\giT,X,)~g{f,X,)\'] < ClE 



Xf — Xf I 



b{X,)ds- 



a{X,)dWs\\ 



where |r - < T|t - f|, E [\\Xf - XfW^] < Cl/i by Proposition 4.1, and 



E 



II / b{X,)ds 



a{X,)dWs\\' 



< E[a|r-f[ 



by Doob's inequality, (HL) and Proposition 4.1 again. 

In order to prove (3.2) of Proposition 3.1, we need the following easy Lemma. 



□ 
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Lemma 4.1. Let (HL) hold. Then, 

max(||Ft.|| + V7^||ZtJ|) <a and \\Y\\s- + \\Z4w + \\Z\\u'2 < Cl 



(4.4) 



Proof. The first bound follows from the same arguments as in the proof of Lemma 3.3 in [4], 
after noticing that the boundedness assumption on h and a can be relaxed for our result. Since, 
by (2.6), 

on combining Jensen's inequality with (HL), the first inequality of (4.4) and Proposi- 

tion 4.1 imply that 

supE[|yi|2] < 2maxE[|Ft,^j2] ^2/i2niaxE[/(Xt,,yt,,Ztj2] < . (4.5) 

i^^j" i<n i<n 

Applying It6's Lemma to F^, using the inequality ab < + b"^ for a, 5 e M, (HL), (4.5) and 
Proposition 4.1 then leads to 



\ZA\'ds 



tAf 



for all a > 0. By Remark 4.2, this shows that 



\ZMs)\\^ds 



E 



\Z, 



4,{s) 



< 



\zJrds 



-1 , „-1ti 



< CL[l + a + a +a 



iZJl'ds 



Thus, taking a large enough, but depending only on L, and recalling Remark 2.2 leads to the 
required bound for II^H-h^ and ||Z0||7^2. The bound on Hi^H^s is then easily deduced from its 
dynamics, Burkholder-Davis-Gundy's inequality, (HL), (4.5) and Proposition 4.1. □ 

Proof of (3.2) of Proposition 3.1. Applying Proposition 4.2 to 9 := t+ A f and recalling 
Remark 2.2 leads to 

Err(/i)2^^- < Cl {h + E [\Yr^^r - n+ArP] + niY)l, + 7^(Z)^.) . 

It remains to show that 









2 




+ E 


lr<rE 


J \\Zs\\ds 1 




) 













(4.6) 
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Since / is L-Lipschitz continuous under (HL), we can find an K^'-valued adapted process x 
which is bounded by L and satisfies 



/(^0(s),^0(s),^</,(s)) = /(-'^0(s),i0(s),O) + (X0(s),^</>(s)) (4.7) 



on [0,T]. Set 



where £ stands for the usual Doleans-Dade exponential martingale, and define Q ~ P by 
dQ/dP = Ht- It follows from Girsanov's theorem that 

Ja 

is a Q-Brownian motion. Now, observe that, by (4.7) and (2.6), 

Yt = g{T,Xr)+ f f{X„Y,,Z,)ds~ [ ZsdWf (4.8) 
Jtr\T Jt/\T 

Yt = g{f,Xr)+ f (f(X^^s),Y^is),Z^is))-lr^<s{XHs),Z,))ds- f ZsdW^ .{4.9) 



t/\f 

In view of (4.7), (4.8), (4.9), it then suffices to show that 

< Cl (h + E 

2" 



E 



E 



[g{T,Xr) ~ g{T,Xr) I J^r+Ar]^ 

/(X0(3),F0(s),O)ds I J^r^ 



[S,L\r-f \ \Tr^r.,\ 



(4.10) 



lr+<rE' 



lr+<rE'' 



< E 



{x<)>{s)iZ^i^s) - Zs)ds I 



{\T-f\+h) I Tr^^f^ , (4.11) 

< Cl/i , (4.12) 



E 



i.<.+eQ 



J{Xs,Ys,Zs)ds I Tf 



< Cl + E 



E[CL|r-f| I Tr^r^f]' 



ClE 



lr<rE 



\Zs\\ds I 



(4.13) 



We start with the first term. By using (HL), applying Ito's Lemma to [g{t,Xt))t>o between t 
and r, using Proposition 4.1, the bound on x well as standard estimates (recall (Hg) and 
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Proposition 4.1), we easily check that on {t+ > r} C {t > t} 



J (1t+<s<t(X0(s)O-*, Ag) + A?) {s,Xs)ds I 



< Cl \\X^-X4+E[^L\T+-f\\Tr] . 

Similarly, on {r+ < f}, 

|eQ [g{T+, Xr^) - g{f, X^) I Tr;\ \ < Cl - 1| + E \r+ - f \ \ Tr^\ . 

We then conclude the proof of (4.10) by appealing to (HL) and Proposition 4.1 to obtain 



E 



I ~ 1 1 2 II - 1 1 2 



-E 



\g{T+,Xr^) - g{T,Xr)\ 



< CLh, 



recall that Q < t+ — t < h. 

The second term (4.11) is controlled by appealing to (HL), Lemma 4.1 and Proposition 4.1, 
recall that t+ — t < h. Concerning the third term (4.12), we observe that {r+ < s} = {t < 
<p{s)} G ^0(s) and that {f > s} = {f > (/)(s)} G ^<p{s)- It then follows from (2.7) that, on 
{r+ < f}, 



:E 



E 



(X0(s) I ^0(s) ~ Zi,)ds I Tr+r\f 

Hs{X<li{s),Z^<^s) - Zs)ds I Tr+ 
f<p{s) + h 



H. 



m \ x,p{s) 



Zudu — Zs ) ds I JF^, 



(i?s — H^j^s)) \X(j>{s) I ^0(s) — Zs) ds I 



and, since r and t+ take values in tt. 



4>{s)+h 



Zudu — Zs ) ds = . 



» 



On the other hand, the Cauchy-Schwartz inequality and the boundedness of x imply that 



E 



{Hs - H^(s)) \ X<P(s) , Z^(s) - Zs) ds I J^r+Ar 



< Cl 



E 


/' 


{Hs — Hff 


^(s))'^ds J-T^/\f 










E 


/' 


ll^</.(s) - 


ZsW'^ds 1 J- r_|_AT 




_Jt+ 







E 
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Recalling Lemma 4.1 and combining the above inequalities leads to (4.12). 

The last term (4.13) is easily controlled by using (HL), Remark 2.2, and Proposition 4.1. □ 

5. Exit time approximation error: Proof of Theorem 3.1 

In this section, we provide the proof of Theorem 3.1. We start with a partial argument which 
essentially allows to reduce to the case where to = 1, i.e. O has no corners, by working separately 
on the exit times of the different domains O^: 

t[ := M{t en : 3s<t s.t. X,, ^ O^} A T and M{t e tt : Xt^O'^} AT . 

We shall prove below the following Proposition. 



Proposition 5.1. Assume that (HL), (Dl) and (C) hold. Then, for each e > 0, 

n2" 



E 



E 



\fl<£< 



(5.1) 



It implies the statements of Theorem 3.1. 

Proof of Theorem 3.1. Since = min^<„i and f = min^<„i f ^, we have 



E 



which combined with (5.1) leads to 



E 



E[|r-f| I .f,+Ar]' 



(5.2) 



since |t_(- — t\ < h. This leads to the second assertion of Theorem 3.1. On the other hand, given 
a positive random variable ^ satisfying E [^p] < C£ for all p > 1, we deduce from Holder's 
inequality that 



E |t - f I I Tr+^,] 



\t - t\ 1-^ I Tr^Af 



2(l-e) 



< a T^' E [|r -f\\T-, 



r+ArJ 



l2(l-e) 



and 



E 



eE[^ |T-f| I Tr^A,]' 



\t-t\ I Tr.Af] 



1-e 



In view of (5.2), this leads to the first assertion of Theorem 3.1, after possibly changing e. □ 

The rest of this section is devoted to the proof of (5.1) for some fixed £. We first provide an 
a-priori control on the difference between and f^. We use the standard idea that consists in 
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introducing a test function on which we can apply ItS's Lemma between and f so that the 
Lebesgue integral term provides an upper bound for the difference between these two times, see 
e.g. Lemma 3.1 Chapter 3 in [9] for an application to the construction of upper bounds for the 
moments of the first exit time of a uniformly elliptic diffusion from a bounded domain. 

To this end, we introduce the family of test functions 

Ft := , 1 < £ < m , 

for some 7 > to be fixed below. Here, di is a C'^iW'-) function which coincides with the algebraic 
distance to dO^ on a neighborhood of dO^ and such that 

O'^ ■.= {x^m.'^ : dt{x)>Q} and dO^ -.^ {x : di{x) ^ Q} . 

The existence of such a map is guaranteed by the smoothness assumption (Dl), see e.g. [11]. 
Observe that, after possibly changing L and considering a suitable extension of d^ outside of a 
neighbourhood of the compact boundary 90^, we can assume that 

\\dt\\ + \\Ddi,\\ + \\D'^dt\\ < L onR'^ . (5.3) 

Observe that 

CFe = l^[{2{b,ne)+T:r[aD^de])de + TT[a{neyne]] (5.4) 

where tii :— Ddg coincides with the unit inward normal for x G dO^, recall (Dl). 

In view of (HL), (Dl), (5.3) and (C), there is some Cl > such that, for each 1 < £ <m, 

£Fi> -{-CLde + nea{ni)*) >1 and ne aim)* >L-^/2 on B{dO\r) (5.5) 

7 

if we choose r > and 7 > small enough, but depending only on L. For later use, also observe 
that, after possibly changing r, one can actually choose it such that 

ni{x) a{y)m{x)* > L-^ /2 for all x, y <a B{dO'^ ,r) s.t. \\x - ij\\ < r . (5.6) 

We now fix r, 7 > such that (5.5) and (5.6) hold and define the sets 

At := {X, eS(aO^r), VsG [f^r^]}, :={|rf,(X,|)|</i^-''} 

At := {X, e B{dO',r) , V .s e [rj,?^]} , Bt {\dt{X^,)\ < /i*""} , 

for some 77 e (0,1/4) to be chosen later on. Observe that At (resp. A() is well defined on 
{f^ < r^} (resp. {r^ < r^}). 

We can now provide our first control on \t^ — f^\. Recall that ^£ (^l if it does not depend on 
some extra parameter e) denotes a positive random variable whose value may change from line 
to line but satisfies E [|^£|p] < C^'^ for all p > 1. 
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E 



>T<!} 



(r-T|)2P {At^Bif I T,. 



1-e 



L{r|<rn 



for each 1 < i < m. 

Proof. 1. We first work on the event {t^ > f^}. It follows from (5.5) and Ito's Lemma that 



E [tI - I J^fi] < E 



lAenBt / CFt{Xs)ds I TfH 



< 



+ iT- f') P [{Ai n Bif I T^,] 



Ia^hbA CFe{Xs)ds+ DFeiX,)a{X,)dwA \ T^i 



lA.nB, / DFeiXsMXs)dW, I 



E 



(T-fOPp.ns,)^ I .F^.] 



< 7" E 



idj{X^c ) - dj{Xfe))lAtnBt I -^T* 



E 



'-{AenBe)- 



DF(iX,)a{X,)dWs I .F^* 



+ iT~f')V[iAenBir 

where, by Holder's and Burkholder-Davis-Gundy's inequality, the Lipschitz continuity of a and 
DFi (see (HL) and (5.3)) and Proposition 4.1, 



< eLiT-f'r^p[iA,nB,r\T,.]' 



MA,nB,y / DFe{Xs)aiX,)dWs \ T^^^ 



for all £ G (0, 1). We now recall that \di(X^i: )| < ''on B^, which implies 



E 



{d\(X^t ) — SfXX^iy)\Air\Bi \ ^fi^ < E d^i{X^t )lAinBi \ 



In view of the above inequalities, this provides the required estimate on the event set {t^ > f^} 
since rj < 1/4. 

2. We now work on the event {r^ < f^}. By Proposition 4.1, 



E 



'-ArnB, 



C^^MFi{X,)-C''^FiiX,) 



ds I 
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with the notation CVFi := dtFi + {b{y), DFi) + ^Tr [a{y)D^Fi] , so that C^^Fe{Xs) = CFi{Xs). 
Arguing as above, it follows that, on {f^ > r^}. 



E 



E 



/ ^ DFt{X,)a{X^^,))dWs I 



(T-t|)P {A,^B,f I 



□ 



It remains to control the different terms that appear in the upper bound of Lemma 5.1. 
For notational convenience, we now introduce the sets (recall that < < 1/4) 

Ei := {deiX^e) < h^"^} and Ee := {di{X^^) < hi"^} , 1 < £ < m . 
Remark 5.1. Observe that 

P [E^e n {?'' < r|}] < P [E'i n {r^ < T}] < P UdiiX^i) - di{X^i) > h^-"} n {f^ < T} 

since de{Xfe) < on {f^ < T}. Using (5.3), Tchebychev's inequality and Proposition 4.1, we 
then deduce that, for each e G (0, 1), there is C£ > such that 

p[s^"n{f^ <<}] < cih^-' . 

Similarly, if denotes the first exit time of {t,Xt)t>Q from [0,T) x O^, we have 



[E^,n{f'>ri}] < 



1 



1. 



{d,ix,) - d,{x,{) > n {d,ix. ) < -h-^-^} n M < t} 



where the last inequality follows from Tchebychev's inequality. Proposition 4.1 and the fact 
that — < h. Note that the term d({X^i^) — di{X^e) could be controlled by Bernstein 
type inequalities in order to avoid the explosion of the constant with e. However, to the best 
of our knowledge, such inequalities are not available in the existing literature for the term 
di{X^i) — di{X^i) and Tchebychev's inequality remains the most natural tool to apply here. 
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Combining the above Remark with the next two technical Lemmas allows to control the right 
hand-side terms in the upper bound of Lemma 5.1. Thus, the statement of Proposition 5.1 is a 
direct consequence of Lemma 5.1 combined with Remark 5.1, Lemma 5.2 and Lemma 5.3 below, 
appHed for 77 small enough. 

Lemma 5.2. Assume that (HL), (Dl) and (C) hold. Then, for each e G (0, 1), 

P [A'i I T,.] lE,n{rl>,^} + P [Ai I Tr^^ lB,n{<<.n < ^ , V £ < m . (5.7) 

Lemma 5.3. Assume that (HL), (Dl) and (C) hold. Then, for each e G (0, 1), 



\Ai n Bl I Tf^\ \Etr\{r'.>f^'} 



A,, n I 



T - A ri 



, V£ < TO 



(5.8) 



Proof of Lemma 5.2. 1. We first prove the bound for the first term. Let V be defined by 
Vt := di(Xfej^.i) for t > and let be the first time when V reaches y G M. Using A'^ = 
n > ^"-j U < d"-}), we deduce that on {t| > f'^} n Ee 



' \A\ I T^t\ < P [?9" > r I Tri] 



{ sup \di{Xs)\ > r} n {t^ < T} \ T^, 



where, by (5.3), Tchebychev's inequality and Proposition 4.1, on {t^ > r^} C {r^ > r^}. 



{ sup \d,{X,)\>r}n{T' <T}\T^. 



< r-^E 

< ah 



sup \di{X,) - diiX,,)\' \ 
s6[t^t|] 



recall that — < h. It remains to provide a suitable bound for P > \ T=ft\ . From now 
on, we assume, without loss of generality, that 



Set ■d:=-d^ K ■&''. Thanks to (C) and (HL), we can define Q ~ P by the density 

H = £r^+.o (^-^ lE^s>f'{ni(7){X,){{ntanl)(X,))-^Cd,{X,)dWs 
Let 

= VF + / {iMay{Xs){{n,anl){Xs))-^Cdi{Xs)ds 



(5.9) 
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be the Brownian motion associated to Q by Girsanov's Theorem. We have 



Set 



At := 



+ / ni{Xs)a{Xs)dW2 on Ee 



r'+t 



By the Dambis-Dubins-Schwarz theorem, see Theorem 4.6 Chapter 3 in [18], there exists a one 
dimensional Q-Brownian motion Z such that 



14a^ = + Za,,, on n {r| > f*^} = {V^ < /i^"" , r| > f'} . 
This imphes that 

Q > -^r^] < on n {r| > f^} , 

see e.g. Exercise 8.13 Chapter 2.8 in [18]. We conclude by using Holder's inequality and (5.3). 
2. The bound for the second term in (5.7) is derived similarly. We now write 

Vt (X,. +J ,t>0. 

As above, we denote by -d^ the first time when V reaches y G M and observe that, by (5.9), 



A'i I 



< 



sup \d,{Xs) - di{Xfi)\ > /i^-" I 



where f ^ := t£ + -d' 



-h2- 



and, by (5.3), Tchebychev's inequality and Proposition 4.1, 



sup \dt{Xs) - di{Xf,)\ > h-^-'^ I 



sG[f',f'+h] 



< ah. 



In order to bound the term 
enough. 



> I T^i 

+ 



, we observe that (5.6) imply that, for h small 



\\m{Xs)a{X^(,))\\> L-^ /V2 on n {5 e [r|, 0^]} n {||X - || < r} , 
where 9^ := mi{t > : Xt ^ B{dO^ ,r)} A T. Moreover, it follows from Proposition 4.1 that 



sup \\Xs - X^(^s)\\ > r 

s<T 



< Cl h . 



Up to obvious modifications, this allows us to reproduce the arguments of Step 1 on the event 
set Ep. □ 
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Proof of Lemma 5.3. We only prove the bound for the first term. The second one can be 
derived from similar arguments (see step 2 in the proof of Lemma 5.2). We use the notations of 
the proof of Lemma 5.2. We first observe that, on Ei n {t^ > t^}, 

F[AenB'g \ Tfi] < ¥[Aer\{'d" > (T-f'^)} \ Tft] 

{t' <T}n sup \diiXs) ~ de{X,e)\ > I T.e 

se[r«,T{] 



< 



Atn{ min Za, > -ft-^"''} I j^^, 
te[o,T— ff] 



ah , 



where the second inequality follows from Tchebychev's inequality, (HL) and Proposition 4.1, 
recall that — < h. Using Holder's inequality, we then observe that 



Ai,r\{ min Za, > -/i^-''} | jT-^ 



Ai,r\{ min Za, > -/i^"''} I jr ^ 



Since, by (5.6), 

At-t^ >{T~ fO(2L)-i on At n > (T - f^} n {f^ < <} C n {f^ < t| = T} , 
we deduce from Chapter 2 of [18] that, on n {f^ < r^}. 



n { min Za^ > -/i^-''} | jr^^ 

te[0,T-r''] 



< 



min Zt > I jF-, 



< Cl (r-f^)-2/i2-'' 



We conclude by combining the above estimates. 



□ 



6. Regularity of the BSDE and the related PDE 

6.1. Interpretation in terms of parabolic semilinear PDEs with 
Dirichlet boundary conditions 

In this section, we denote by X*'^ the solution of (1.1) with initial condition x <E O ai time t <T. 
We also denote by r*'^ the first exit time of (s,X*'^)s>t from O x [0,r) and write (y*'^,Z*'^) 
for the solution of (1.2) with (X*'^,t*'^) in place of {X,t). 

As usual the deterministic function {t, x) ^ D ^ u(t, x) := Y^'^ can be related to the semiHnear 
parabolic equation 

= -Cu{t, x) - fix, uit, x), Du{t, x)a{x)) , {t,x)eOx [0, T), 

I (6-1) 
u\dpD = g ■ 
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where we recall that C denotes the Dynkin operator associated to the diffusion X, £V := 
dtijj + {b,D^) + iTr [aD^^jj] with a := aa* , and dpD := ([0,T) x dO) U ({T} x O) is the 
parabolic boundary of D. 

Proposition 6.1. Let (HL), (Dl), (D2), (C) and (Hg) hold. Then the function u has 
linear growth and is the unique continuous viscosity solution of (6.1) in the class of continuous 
solutions with polynomial growth. 

A similar result is proved in [6] but in the elliptic case. For the sake of completeness, we provide 
a slightly different complete proof of the viscosity property in the Appendix, where the standard 
associated comparison result leading to uniqueness is also stated. 



6.2. Boundary modulus of continuity 

Adapting some barrier techniques for PDEs, we first prove the following bound for the modulus 
of continuity on the boundary. 

Proposition 6.2. Let (HL), (Dl), (D2), (C) and (Hg) hold. Then, there is Cl > such 
that for all {tQ,xo) £ [0,T) x dO, 

lim K^o.2/)-"(;o^^o)l ^c,. (6.2) 

In particular, if the gradient of u exists at {to,XQ), it is uniformly bounded. 

Proof Let {to, xq) e [0,T)xdO and A -.^ [to,T)xJ\f, where A/" C O is an open set and xq e dJ\f. 
We only show that, for all y £ M, 

uito,y)-u{t,,xo) ^ 

\\y - Ml 

The lower bound is obtained similarly. By (D2), there is e > and a family (ei)ig|i d| such that 
xq + eei € N for all i £ and span(e,;, i £ |l,(i]) = M.'^. Thus, (6.2) implies the statement 

concerning the gradient, whenever it is well defined. We now prove (6.3). 

1. Assume that there exists a smooth function : A M. with first derivative bounded by Cl 
such that 

(a) ip>u on dpA := {[to, T) x dU) U ({T} x Sf). 

(b) Cil;{t, x) + f{x, ip{t, x), Dip{t, x)a{x)) < for {t, x) £ A. 

(c) ilj{to,xo) ^ u{to,xo) = g{to,xo). 
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Using Proposition 6.1 and a standard maximum principle, see Lemma A. 2 in the Appendix, we 
then derive that u < %jj on A. In view of (c) this yields 

n n < n n <Cl , Vy e TV \ {xo\ • 

lly-2;o|| \\y-xo\\ 

2. It remains to construct a smooth function satisfying (a), (b) and (c). Recall that the spatial 
boundary dO is compact. Since u is continuous on D, see Proposition 6.1, the compactness 
assumption (Dl) ensures the uniform boundedness of u in a neighborhood of [0,T] x dO. 
We specify the construction of the barrier function only for xq G dO\B{C,L~^). Indeed, for 
xo G B{C,L~^), assumption (C) ensures that the diffusion coefficient is uniformly eUiptic in a 
neighborhood oi xq. The expression of the barriers below can then be simplified. Namely, we do 
not need the additional localization with the cone, i.e. we can take k = in (6.6) below. 

Let y := y{xf)) he the point ofO'^ associated to xq by the exterior sphere property, see (D2). Set 
r := r{xo) ~ \\y{xo) — a;o||. Recall that, by assumption, B B{y, r) satisfies B O O — {xq} . 
It follows from (HL) and (C) that 

{a{x)n{xo),n{xo)) > on the set Vi -.^ {x e O : \\x ~ xo\\ < Vl} (6.4) 

for some rjL > small enough, but depending only on L. 
For X € O, we now set 

dB{x) := d{x, dB) = \\x — y\\ — r 

so that e C2((5) with 

( \ n2^ ( \ {x-y)*{x-y) 
DdB(x) = -j^ 7- , D dB(x) = -j^ 7 (6.5) 

\\x-y\\ \\x-y\\ \\x-yr 

where Id denotes the identity matrix of M'' . We now introduce a cone 

K := {x e R"^ : (a; - y, n{xo)) > cos(6')||a; - 6 e [0, tt/2] 

and 

r>2 := {x e O : dsix) < 6} , 6>0, 

where S < Sl small enough to ensure I?2 C We finally set A/" := O n /C n I?2 and define the 
barrier function by 

^(t, x) := 5(t, x) + Aaiifiix)'/' - S'/^) + k{x - y, n(xo)) (^1 - ^^^^^^^f^) (6-6) 

for (t,x) e [to,T] X JV, where ip{x) := S + dsix) . for some (a, k) G (0,oo)^ to be chosen later 
on. 
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ao 



K. n(x\ 



Figure 1. Domain for the barrier 



2.b. Since xq — y <E span(n(xo)), il}{tQ,XQ) — u{tQ,xo) — g{tQ,xo), so that (c) is satisfied. 
2.C. Recall from the beginning of Step 2. that 

M := sup |M(t,a;)|V sup \g{t,x)\ <oo. (6.7) 

{t,x)e[tQ,T]xt>i {t,x)e[to,T]xVi 

On dOndM, Tp{t,x) > g{t,x). On dV2ndAf, tp{t,x) > -M + Aa{2^/^ - 1)5^/^ . Thus, for 

M , , 

" ^ 2(21/2-1)51/2' (6.8) 

one has ip{t, x) > u{t, x) for {t, x) G [to, T] x 9X>2 n dAf. 
On a/C n 97V, we have 

ipit,x) > -M + Kcos{0)\\x ~ y\\{l - coa{0)) > -M + nr cos{6){l - cos{6)) . 

Hence, for 

2M 

rcos(0)(l-cos(0))' (^'^^ 

we obtain that ^p(t, x) > u(t, x) V(i, x) G [to, T] x dIC n OA/". This concludes the proof of (a). 
2.d. It remains to show that satisfies (b). Set 



V \\x~y\\ J 



r( 

and observe that, for some C < Cl, 

|lOT(x)|j<C, \\D'r{x)\\<C/r (6.10) 
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uniformly in x G Af. Define, 

e(t,x) := Cil;{t,x) + f{x,ilj{t,x),Dij{t,x)a{x)) 

< c(i + M + + „(i + .-.))- 1 fwiif^. ^) ,i.r"' 

r 

< C{l + M + .{1 + r-i)) - f ^(x)-^/^ - C{1 + r-^Mx)^ 

recall (Hg), (6.5), (6.7) and (6.10). For a suitable angle of the cone 9, we shall show below that 
we can find C > such that C^^ < Cl and 

a(x)^^, „^~^,A >C , Vx G . (6.11) 



lla^ - y\\ ' 11-^; - y\ 

Recalling that </?(a;) < 26 for x G A/" C I?2 , we get 

e(t, x) < C(A/ + >i{l + r-^))- ^ip{x)-^^^ (C - 2C(1 + r^i);^) . 

ForfS := {l/A)C{C{l+r-^))-^ A5l > 0, we then havee(t,a;) < C{M + K{l + r-^))-Ca2-i . 
It is then clear that (a, k) can be chosen in order to satisfy (6.8), (6.9) and so that 0(t, x) < 0. 
This shows (b). 

It remains to prove (6.11). This is done by suitably choosing the angle of the cone K.. Let 
Z G 3-6(0,1) be such that Z + y E IC. Introduce the basis (n(a;o), (?T.^(xo))ieji^c;-i]) where 
(n^(xo))ig|i_d-i] is an orthonormal basis of {n{xo)}-^ for the euclidean scalar product. Let 
(A)ie[o,ci-i] denote the coefficients of Z in this basis, i.e. Z = f3on{xo) + X^fji Pinj^ixo) . One 
has, for all x & J\f, 

d-l 

{a{x)Z,Z) = /3o(a(.T)n(a;o),n(xo)) + 2^/3o/3i(a(a;)n(xo),n^(xo)) 

i=l 

d-l d-l 

+ {a{x)^(3,nf{xo),^f3,ni{xo)) 

i=l i=l 

d-l 

> PQ{a{x)n{xo),n{xQ)) + 2^ f3af3i{a{x)n{xa),nj'{xo)) . 

1=1 

Since Z + y E K. and \\Z\\ = 1, we must have /3o > cos0, by definition of IC, and therefore 
< sin(6') for alH G |1, d — 1]. Hence, (6.4) and the above equation leads to 

T-l 

{a{x)Z,Z) > cos^{0)— 2(d- l)sin(6') sup ||a(a;)|| , Vx G TV. 

2 xeN' 

This yields (6.11) with C = - — for 9 small enough. □ 
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6.3. Representation and weak regularity of the gradient in the regular 
uniformly elliptic case 

In the section, we strengthen the initial assumptions and work under: 

(D'): O is a bounded domain satisfying (Dl) and (D2) for the constant L. 

(C): a is uniformly elliptic with ellipticity constant L~^. 

(H'): the coefficients b, a, f and g satisfy (Hg)-(HL) and are uniformly C'^{D). 

Prom now on, given a matrix M , we denote by M'^ its j-th column, viewed as a column vector. 

Proposition 6.3 (Representation of the gradient). Let the conditions (D'), (C) and (H') 
hold. Then, u e C^{D) n C^^^{D), Du e C^{D) and for all (t, x) <E D 



Du{t,x) = E 



(6.12) 



where V X^-^ is the first variation process of X^-^ : 

VX*'^ ^Id + Yl / Dcr-'{Xl-'=)\/Xl^'=dW.i + / L'6(X*'^)VX*'^dw ,s>t, 

and y*'^ is defined by 

V*'- := cxp (^^' dyfielndv + a./(e*--)dt¥„ - \\dJ{QYWdv^ ,s>t, 

with e*'"" = (x*^^,y*'^,z*'^). 

Proof. The result is obvious for (t, x) e dD. We then assume from now on that (t, x) G D. We 
derive from Theorems 12.16 and 12.10 in [21] and the definition of Holder spaces at p. 46 of this 
reference that Du G C^{D). Let us consider the systems of differential equations obtained by 
formally differentiating the PDE (6.1) w.r.t. (x')igji_(;]- For i = 1, . . . , c?, we have 

= dtv' + {b + (7*DJ{Q) + ]^D^^a■\Dv') + ]ph:[aD'^v'] (6.13) 



{D^.b' + Dyf{Q) + {DJ{Q), D,.a')) + D,.f{Q) + ^ h^''^ , 

d 

where /i'''' ^ {D^.b'' + {DJ{e),D^^a^'')) D^kU + ^ D^.a''^ D^k^iu 

1=1 



and Q{t,x) = {x,u{t,x), Dua{t,x)). 

Given n large enough, set £>„ := {x + B{0,n-^), x £ O'Y C O, Tn := T - n-^ > and 
Dn ■= [0,Tn) X On- Note that by construction On satisfies a uniform exterior sphere property 
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(with radius l/2n). Then, the PDE (6.13) on D„ with the boundary condition D^iu on dpDn = 
([0,r„)xae'„)U({r„}xO„) admits a unique C°{D„)r\C^'^{Dn) solution <, see Theorem 12.22 
in [21]. Using the maximum principle, we can then identify D^iU and on l)„ by considering 
the PDE satisfied by e~^{u{-,x+eei) — u{-, x)) — vl^(-,x) on Z)„. Here, is the «-th canonical basis 
vector of M'', see e.g. Theorem 10 Chapter 3 in [10]. In particular, Du G C°{Dn)r]C^''^{Dn). By a 
usual localization argument, we then deduce from ItS's Lemma applied to Du{-, X*'^)VX*'^F*'^, 
with {t,x) G Dn, that 

Du{t, x) = E Du{Tn,Xl^)VX'^'^V^f + " (9^/ (6*'^) VX*'^V;*'^ds 

where r„ :— inf{s G [t,T„] : (s,X*'^) ^ Dn}- Observe that lim„ r„ = t P — a.s. by continuity 
of X. We then derive the statement of the Proposition by sending n — > cx), using the a-priori 
smoothness of u, Du e C"{D), and the dominated convergence theorem. □ 

Remark 6.1. Note that the various localizations in the previous proof are needed because we 
do not assume any compatibility condition on the parabolic boundary, i.e. Cg + f{-, g, <jDg) = 
on dpD. Otherwise, Theorem 12.14 in [21] would give u € C^''^[D) which would allow to avoid 
the introduction of the subdomains On- 

Observe that, by Proposition 6.2 and the continuity of Du stated in Proposition 6.3, we have 
||£>w(t*'^, X'i'^^)!! < Cl- The representation (6.12) and standard estimates then give ||£'u||oo,i5 — 
Cl- 

Corollary 6.1. Let (D'), (C) arid (H') hold. Then, \\Du\\^ q < Cl- 

We can now prove Theorem 3.2 under the conditions (D'), (C) and (H'). 

Corollary 6.2. Theorem 3.2 holds under the conditions (D'), (C) and (H'). 

Proof. 1. Proof of (3.4) and (3.5). Recalling that u G C'^^'^{D) n C^iD), see Proposition 
6.3, we deduce from a standard verification argument that Z ~ Du{-,X)a{X). Set {\7X,V) 
(yXO.Xo^yO^Xo-^ and observe that (VX*-'^*, K*'^*) = (V X^V X^r\VsVt-^) for s>t,hy the flow 
property. Thus, by Proposition 6.3, 



E 



Du(t, XrjWXrVr + d,.f (6,) VX.Vsds \ Tt 



G{Xt){VXtYt)-^ ,t<T. (6.14) 



It then follows from Proposition 6.2 (boundedness of the gradient of u), (HL) and stan- 
dard estimates that supt<^ \\Zt\\ < £,L. This readily implies (3.5), i.e. E // \\Zs\\Pds | J^e < 
E[^£|i9 — 0| I Tg], p = 1,2. By Burkholder-Davis-Gundy's inequality, (HL) and Proposition 
4.1, this also yields E [sup^gfe^^j IF* - YeM < E [^l \d - 6*1^], p > 1 
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2. Proof of (3.6). By the same arguments as above, we first obtain that \u{t,x) — u{t,x')\ < 
Cl\x - x'\. Moreover, for t<t' <T, 

u(t, x) — u{t' , x) ~ Y^'^ — u{t\ x) = y/'^ — Yj*'^ + u{t' , Xli^) — u{t\ x) . 

The Lipschitz continuity of u in space (Corollary 6.1) and standard estimates on SDEs imply 
that \E[u{t',Xl:''y u{t',x)]\ < CL\t-t'\i. On the other hand, E [1^/'"= - <CL(t'-t), 
by the above estimate. 

3. Proof of (3.3). The bound on TZ{Y)^2 follows from (3.4). Using (6.14) and exactly the same 
arguments as in the proof of Proposition 4.5 in [3], see also [2.3], we deduce that 



\Zt-ZtM^dt 



<CLh , 



which impHes J27=o ^ \\Zt - ZtM^dt 



< Cl h since Z is the best approximation of Z in 



L^{fl X [0,T]) by an element ofH^ which is constant on each time interval [ti,i,;+i). 



□ 



6.4. Regularization procedure: proof of Theorem 3.2 in the general 
case 



Step 1. Truncation of the domain: We first prove that Theorem 3.2 holds under the condi- 
tions (Dl), (D2), (C) and (H'). 

Let (j) he a density function with compact support on R'^. Given e > 0, we define := 
£~'^4>{£~^ ■)-k{dAdi,-i)~^ where d^^i denotes the algebraic distance to dB{XQ, e~^) and ★ denotes 
the convolution. Set := {x e R'^ : Ae(x) > 0} and := [0,T) x C^- It follows from the 
compact boundary assumption that dO C Oe, for e small enough. Note that Oe is bounded, 
even if O is not. Let (V^ , Z^) be defined as in (1.2) with Oe in place of O and be the first 
exit time of (-jX) from De- Observe that, by continuity of X, r"^ t P — a.s. Since, by (Hg), 
(HL) and Theorem 1.5 in [26], 



\Y -Y' 



2 

I 52 



\Z-Z^ 



< ClE 



\' + \YJ' 



f{Xs,Ys,Zsfds 



'-)ds 



we deduce from Proposition 4.1 and a dominated convergence argument that — y^"^ 



52 ' 



\Z- 



0. Since the domain satisfies (D'), we can apply Corollary 6.2 to {Y^ , Z^). Recalling 
that the associated constants depend only on L and are uniform in e, we thus obtain the required 
controls on {Y,Z). Let be the solution of (6.1) associated to Dg. The above stability result, 
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applied to general initial conditions, implies that u'^ ^ u pointwise on D. Corollary 6.2 thus 
implies that u satisfies (3.6). 

Step 2. Regularization of the coefficients: We now prove that Theorem 3.2 holds under 
the conditions (Dl), (D2), (C), (HL) and (Hg). 

For e > 0, define b^, and by 



{be,(yeJe){x,V,z) := {b, a, f) e' 



'2d+l, 



where is a C°° density function with compact support on 



FBSDE 



Let us consider the 



(6.15) 



where {Wt)t>o is an additional d-dimensional Brownian motion independent of W and 

:= M{s > : (s,Xf) ^ D} . 

This system satisfies the conditions of Step 1. Therefore, the estimates of Theorem 3.2 can be 
applied to (V^, Z^). Note that the associated constant depends only on L and are uniform in e. 
Moreover, it follows from (HL) and Theorem 1.5 in ['2f>] that 



\Y-Y'\\i 



52 



\Z-Z'f^2 < C'lE 



|g(r,X.)-5(T^X,-.)r+ \\X,~Xirds 



{\f{X,,,Y,,Zs)\ + \.UXI,Y,,Zs)\)^ds 



L e . 



Clearly, X^ X inS^. Since / and g are Lipschitz continuous, / and have linear growth and 
(X, X^, Y, Z) is bounded in S'^ x S"^ x S'^ x Ti'^ , it suffices to check that ^ r in probability 
to obtain the required controls on {Y,Z). This is implied by the non-characteristic boundary 
condition of (C), see e.g. the proof of Proposition 3 in [17]. The control (3.6) is obtained by 
arguing as above. □ 

Appendix: Proof of Proposition 6.1 

In the following, we use the notations 

u*{t,x) ~ limsup u(s,y) , u*(i,a:) = liminf u(s,y) , (t, x) € D. 

{s,y)eD~*{t,x) {s,y)eD^(t,x) 



The statement of Proposition 6.1 is a direct consequence of Lemmas A.l and A. 2 below. 
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Lemma A.l. Let the conditions of Proposition 6.1 hold. Then, the function u has linear 
growth andu* (resp. u^,) is a viscosity suhsolution (resp. supersolution) of {6.1) with the terminal 
conditions u* < g (resp. > g) on dpD. 

Proof. 1. The linear growth property property is an immediate consequence of Proposition 4.1. 
2. It remains to prove that u* and u» are respectively sub- and supersolution of (6.1) with the 
boundary conditions u* < g and u* > g on dpD. We concentrate on the supersolution property, 
the subsolution property would be derived similarly. The proof is standard, as usual we argue by 
contradiction. Let {to,XQ) G [0,T] x O and tp he such that = min^j j.-,g£,(u* — (p){t,x) = 
(u* — (p){to, Xq) where the minimum is assumed, w.l.o.g., to be strict on D. Assume that 

{-'C(p{to,XQ) - f{xo, (p{to, Xq), DLp<7{to, xo))) l(^t„^^„)eD + -5)(io,a;o)l(to,xo)6apn ~- < ■ 

Recall from (D2) that if xq G dO then we can find an open ball Bq C such that Bq CiO ^ 
{xo}. If Xq G do, we denote by the algebraic distance to Bq. On D, we set 

d{x) 



(p{t, x) = Lp{t, x) - (Vr - t)lta=T - d{x) ( 1 — ) lxoeaO\B(C,L-i) 



. / N A dpo {x) \ 

-aBo[X) I i I lxo€dOnB{C,L-i} , 



for some 77 > . Observe that {to ,xo) is still a strict minimum of (u* — (^) on 1^, n l) for some open 
neighborhood of {to, xq) on which (dso V d) < rj/2 if xq G dO. Without loss of generality, we 
can then assume that 

u>u^>(p + C on dVy^XD" , (A.16) 

while 

<^<V<5-C on V^C^ dpD , if (to, a;o) € dpD . (A.17) 

Moreover, observe that for F equal to d or dso, D{F{l-F/'q)) = DF{l-2ri-^F) and D'^{F{1- 
F/rf)) = (1 - 2r]-^F)D^F - 2r]-^DF*DF where \\DF\\ = 1. Thus, (C) implies that, for 77 and 
Vrj small enough, 

- C^- /(•, (f, D^a) < -C < on K, n 5 . (A.18) 

Let {tn,Xn)n be a sequence in D n such that (t„, a;„, M(t„, {tQ,XQ,u^,{tQ,XQ)). Let 

(A", y", Z") be the solution of (1.1)-(1.2) associated to the initial conditions (<„, Xn) and define 
On as the first exit time of n VJ, by (•, A"). By applying Ito's Lemma on (p and using (A.17), 
(A.18), (A.16) and the identity u = g on dpD, we get 



^{tn,xn) = -x + u{en,Xl) + {f {Xl\ ^{s , XI') , Dipa{s , XI')) ~ r^,)ds 
D^a{s,X^)dW, , 
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where x is a bounded random variable satisfying x > C P — a.s. and 77 is an adapted process in 
such that 77 > C dt X dP-a.e. Following the standard argument of the proof of Theorem 1.6 in [26], 
we deduce that (p(tn, Xn) < Y^^'^" — C,e~^^ = u{tm Xn) — C,e~^'^ . Since (p{tn, Xn) — u(tn, x„) — > 0, 
this leads to a contradiction. □ 

We now state a comparison theorem for the PDE (6.1). The proof is quite standard, see e.g. [5], 
but we give it for the sake of completeness. 

Lemma A. 2. Let the conditions of Proposition 6.1 hold. Fix to G [0,r) and N C O an 
open set. Let U (resp. V ) he an upper-semicontinuous suhsolution (resp. lower- semicontinuous 
supersolution) with polynomial growth of (6.1) on A := [tQ,T) x H such that V > U on dpA := 
i[to,T) X dAf) U ({T} X JV). Then, V>U onA. 

Proof. Fix p > and observe U and V defined by U{t, x) = U{t, x)eP* and V{t, x) = V{t, x)eP* 
are sub- and supersolution of 

pij{t,x) ~ Cip{t,x) ~ f{x,e-P*^{t,x),e~P*Di;{t,x)a{x)) , (i, x) € [io, T) x .(A.19) 

As usual we argue by contradiction and assume that supjj ,j-)g^(C/(t, a;) — V{t,x)) > 0. Define 

P{t,x) :=e-"*(l + ||a;fP) , {t,x) e A 

for p eW such that {\U{t,x)\ + \V{t,x)\)/{l + ||a;||'') is bounded on A, and k > to be chosen 
later on. For all e > small enough, we can then find {t^, x^) G A such that 

sup {U{t, x) - V{t, x) - 2s/3{t, x)) =: {U{t„ Xe) ~ T>(<e, Xe) - 2eP(te, Xe)) > . (A.20) 
{t.x)eA 

Clearly, (t^^Xe) ^ dpA since U < V on dpA. For n G N*, let {tn,Xn,yn) G [to,T] x TV"^ be a 
maximum point of 

(t/(t, x) - V{t, v) - eiPit, x) + Pit, y)) ~ {\t ~ t,\^ + \\x - x,\\^ + n\\x - yf) . 
It is easy to check, see e.g. Proposition 3.7 in [5], that 

U{t^,Xn)-V{tn,yn) ^ {U~V)%,x^) and |i„-teP + lk„-a;ef + n||x„-y„||2 ^ . (A.21) 

Since (te,Xe) G A, we can assume that (t„,a;„) G A for all n G N*, after possibly passing to 
a subsequence. It then follows from Ishii's Lemma, Theorem 8.3 in [-5], that we can find real 
coefficients a„, 6„ and symmetric matrices Xn and such that 

(a„,p„,A'„) G V^U{tn,Xn) and (&«,?«, 3^n) G 'P^V(tn,yn) , 
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see [5] for the standard notations and P^, where 

p„ 2n(x-„ - tin) + A{xn - a:^£)||a;„ - Xe\\^ + eD/3{tn,Xn) , Qn := 2ri(x„ - y„) - eDj3{tn,yn) 
and 

a„-&„ = 2(i„-t,)+e(at/3(t„,x„) + 9t/3(i„,y„)) , I ° \ < An+n-^^^A^f (A.22) 

with 

( 4/d||.T„ - XeW^ + %[Xn - X,)*{Xn - X^) \ 

where /d is the identity matrix of M''. Since U and V are sub- and supersolution of (A. 19), it 
follows that 

p (u{tn,Xn) - V{tn,yn)^ < 0,^ - bn + {b{Xn),Pn} " (^(yn),?™) + ^Tr [a{Xn)Xn ~ a(yn)>'n] 

+ {f{Xn,U{tn,Xn),e~''*"PnCr{Xn)) ~ /(^n , ^(tn ,?/„), e"''*" 5„Cr(y„ ))) e''* 

We then deduce from (HL), (A.22), (A. 21), and standard computations that 

) — V{tn, y„)^ < L (jj{tn, Xn) — V (trn J/n)^ + 2£ (£/3(te, X^) + L 1 1 ^(Ze )_D/3|| ) + 0„(1) . 

Taking p > 2L and k large enough so that £/3 + i||(T-D/3j| < — ■|cxp(— kT) on A, which is 
possible thanks to (HL), we finally obtain 

]^p{u{tn,Xn) -V{tn,Vn)^ < -K CXp(- KT)e + 0„ (1) , 

which contradicts (A. 20) for n large enough, recall (A. 21). 

□ 
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